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SUMMARY 

A detailed description is given of a method oj 
approximating solutions to the integral equation that 
relates oscillatory or steady lift and downwash dis- 
tributions on finite wings in subsonic flow. The 
method of solution is applicable to general plan 
forms with either curved or straight leading and 
trailing edges. Moreover , it is directly applicable 
to control surfaces such as all-movable tails but 
modifications are needed to apply it to controls in 
general. Applications of the method involve evalua- 
tions of numerous integrals that must be handled by 
numerical procedures but systematic schemes of 
evaluations have been adopted that are well suited 
to the routines of automatic digital computing ma- 
chines. These schemes of evaluation have been in- 
corporated in a program for an IBM 704 electronic 
data processing machine. With this machine , a 
pressure distribution together with such quantities as 
section or total lift and moment coefficients or general- 
ized forces can be determined for a given value of 
frequency arid Mach number and for several (four or 
five) modes of oscillation in about 4 minutes of ma- 
chine time. In the case of steady downwash condi- 
tions corresponding quantities can be obtained in 
about 2 minutes of machine time. 

In order to illustrate applications of the method \ 
results of several calculations are presented . In 
these illustrations total forces and moments are com- 
pared ( 1 ) with results of analytic procedures for a 
circular plan form with steady downwash conditions, 
(2) with results of other theories and with experi- 
ment for a rectangular plan form of aspect ratio l at a 
uniform angle of attack, and (3) with some experi- 
mental results for a rectangular plan form of aspect 
ratio 2 undergoing pitching and flapping oscilla- 
tions. Also included in the illustrations are results 


of flutter calculations compared with experimental 
results for an all-movable control surface of aspect 
ratio 3.50 and for a cantilevered rectangular plan 
form of aspect ratio 5.04. 

INTRODUCTION 

The equations that provide a basis for most of 
the existing theoretical aerodynamics of lifting 
surfaces result from a linearization of boundary- 
value problems for velocity potentials. The 
linearizations, of course, involve restrictions in 
addition to that of neglecting viscosity, which is 
required in formulating the potential equations in 
the first place, but these restrictions are generally 
necessary in order to simplify the problems to 
tractable forms. A main advantage gained through 
linearization is the principle of superposition, 
inherent to linear differential equations, which 
means not only that a given problem may be 
broken down into a category of different problems, 
which separate such parameters as thickness, 
camber, and angle of attack, but, more signifi- 
cantly for present purposes, that the linearized 
boundary-value problem can be converted to an 
integral equation; furthermore, the linearized 
differential equation is satisfied not only by the 
velocity potential but also by the pressure or a 
pressure potential so that the integral equation 
can be formulated on the basis of either a velocity 
potential or a pressure potential. 

Even with the simplification and advantages 
gained through linearization the lifting-surface 
problem for oscillating wings lias been explicitly 
solved for only a very few special wing plan forms. 
These are the wing of infinite aspect ratio (refs. 1 to 
3), the wing for which the flow normal to all edges 
is supersonic (ref. 4), the wing of vanishingly small 

l 
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aspect ratio (refs. 5 and 6), and, for an incompres- 
sible fluid, the wing of circular plan form (refs. 7 to 
9). Not only do these few solutions constitute 
a small category but the prospects of obtaining 
explicit solutions for more general or .arbitrary 
plan forms are extremely dim. In view of com- 
puting equipment that is now available, however, 
there exists an excellent possibility of obtaining 
solutions appropriate to any plan form by approxi- 
mate or numerical procedures that are satis- 
factory for engineering purposes. A reliable 
concept on which to base these procedures under- 
lies such developments as those of Falkner (ref. 10) 
and Multhopp (ref. 11) relating to steady finite 
wings and that of Possio (ref. 12) relating to 
oscillating wings of infinite aspect ratio. 

This concept pertains to the bound ary- value 
problem for the pressure potential when expressed 
as an integral equation that relates a prescribed 
downwash distribution to an unknown lift distri- 
bution. It is simply that from a knowledge of 
the few known explicit solutions to the lifting- 
surface problems the general character of the lift 
distribution for various conditions can be sur- 
mised; hence, by assuming that its general char- 
acter is known, the unknown lift distribution can 
be replaced by a sum of appropriately chosen 
modes of lift functions, each mode being weighted 
by a constant coefficient that must yet be deter- 
mined. Employing this concept or procedure 
then allows the known downwash distribution to 
be expressed as a sum of definite integrals with the 
unknown coefficients appearing as factors of the 
integrals. The definite integrals, however, arc 
extremely unwieldy and, in general, can he evalu- 
ated only by approximate or numerical methods; 
it is in this regard that various lifting-surface 
methods, based on the concept under discussion, 
have their basic differences (i.e., in regard to the 
scheme or method by which integrations arc 
numerically accomplished). Falkner’s procedure, 
which has been extended to an oscillating funic 
wing in an incompressible medium in reference 
13 and to an oscillating finite wing in a compress- 
ible medium in reference 14, involves the use of 
certain numerical integrating schemes which, 
although relatively simple to employ, depend in 
part on two-dimensional aerodynamic considera- 
tions and are difficult to assess with regard to 
accuracy. In Multliopp’s procedure, which lias 
been considered for the oscillating case in refer- 


ences 15 to 19 and others, more straightforward 
integrating schemes are used than are employed 
in Falkner’s method but Multhopp’s method, as 
he implies, is devised on the basis of a compromise 
between accuracy of results and computing effort 
required. It must be recalled, however, that both 
the Falkner and Multhopp methods were formu- 
lated at a time when high-speed computing equip- 
ment was in an embryonic state of development. 
Fortunately, with present-day computing equip- 
ment, it is seldom necessary to sacrifice accuracy 
of results for savings in labor of calculation. 

In the present report a method is described in 
which straightforward integrating schemes that 
arc well suited to the routines of high-speed 
digital computing machines are employed. The 
integrations can be made as accurate as desired; 
the accuracy can be easily assessed by simply 
changing the density of steps. The method 
applies to wings with oscillatory downwash condi- 
tions as well as to wings with steady-state down- 
wash conditions. It also applies, with appropriate 
modifications in chosen modes of lift and regions 
of integration, to supersonic speeds as well as to 
subsonic speeds; however, only subsonic speeds 
arc considered heroin. The procedure is applicable 
to general plan forms with either curved or straight 
leading and trailing edge's. Although it is readily 
applicable to control surfaces such as all-movable 
tails, modifications would be needed to apply it 
to controls in general. 

A brief discussion is presented of the integral 
equation under study; the terms of the equation 
arc examined and put in convenient forms, and 
a suitable working form of the equation is obtained. 
The numerical processes involved in handling 
this working form of the integral equation arc 
discussed in detail. Still left untreated are some 
problems concerning the optimum number and 
location of control points and the degree of con- 
vergence of the process. Results of applying the 
procedure are given for a variety of plan forms 
and conditions. These include determinations of 
aerodynamic forces and moments for some special 
conditions and applications to some flutter 
problems. 

SYMBOLS 

A aspect ratio, 4 P/S 

A n (rj) span wise loading function in lift dis- 

tribution (eq. (13)) 
a velocity of sound, ft /sec 

a nm weighting factor in A n (ij) 
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b local wing scmichord, ft 

bo wing root semichord, ft 

Cl, Qp C m , Qj complex total lift and pitching- 

moment coefficients associated 
with jih mode of oscillation; for 
example, C L ,„= (C L ., J ) r +i(C L ,, J ) ( 
\^L.<,j\,\C M , tj \ magnitudes of complex lift and 
pitching-moment coefficients; for 
example, 

F{y) complete chord wise integral in work- 

ing form of integral equation (eq. 

( 28 )). 

I 4 v complete chordwise integral at sta- 

tion p (eq. (35)) 

fn(rj) chordwise integral which involves 

ft th chordwise term in lift distribu- 
tion (eq. (29)) 

n(x,y,t) vertical displacement of wing surface, 

positive down, ft 

hj(x f y) shape function of vertical displace- 

ment in jth mode 

Ii(k$\y—y\) modified Bessel function of first kind 
of first order 

portions of spanwise integral (eq. 
(32)) 

Ki(kx\y — r}\) modified Bessel function of second 
kind of first order 

K[M,k,x 0 , kernel function of integral equation, 

_j(y—y)] ft " 2 

KlMjkjXo, modified kernel function (eq. (16)) 

Ke(M,k f 9,Tj) modified kernel function in terms of 

angular chordwise coordinate 6 
(see eq. (24)) 

k reduced frequency, 6 0 «/F 

L total lift, positive down, lb 

Li(kx\y—r}\) modified Struve function of first order 

L(6,r}) nondimensional lift function (cq. 

a° }) 

l wing semispan, ft 

K{9) ftth chordwise pressure mode ex- 

pressed in terms of angular chord- 
wise variable 0 (eq. (12)) 

ln(k) nth chordwise pressure mode ex- 

pressed in algebraic form (eq. (14)) 
M Mach number 
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M a 

tolal pitching moment, positive nose 


up, ft-lb 


M«=p ^S&oS|jGr. fi 


local pressure difference between top 

*p(8,n,t) 

and bottom surfaces of wing, posi- 
tive down, lb/sq ft 

<h(t) 

generalized coordinate of jih degree 


of freedom (as a subscript, re- 

placed by a denotes pitching, q j 
replaced by <p denot es flapping, and 
qj replaced by x 2 denotes parabolic 
camber deformation), ft 

S 

area of wing surface, sq ft 

s 

ratio of wing semispan to root semi- 


chord, lfb Q 

t 

time, sec 

V 

velocity of airstream, ft/sec 

w{x,y) 

amplitude function of prescribed 


downwash, w(x,y,t)=w(x,y)e i:at , ft/ 


sec 

x ,i 

x 0 =x—£ 

chordwise coordinates referred to b 0 

y,n 

spanwise coordinates referred to l 

Z 

vertical coordinate (sketch 1) 

/3=-Vl - AP 

5 

spanwise interval between integra- 


tion stations in region of span 
which contains control point 

f 

semispan of region which contains 


control point (sketch 4) 

8 

angular chordwise variable (eq. (7) 


and sketch 2) 

7 

m 

O 

w* 

II 

(£>»— .r)J 

A 

angle of sweep of quarter-chord line, 


measured from perpendicular to 
stream direction, positive clock- 
wise, deg 

\ 

, Tip chord 

lapcr ratio, Root o|iord 


model mass ratio (mass of model 


divided by mass of air in a circular 
cone with local diameter equal to 
local cl lord of model) 


coordinates of leading edge, trailing 


edge, and midchord line, respec- 
t ively 
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P 

<#> 


<f>L .q .Q j 


CO 

0>j 

Indices: 

i,r 

3 

n,m 
x 2 

Of 

<t> 


air density, slugs/cu ft 
generalized coordinate of angular dis- 
placement in flapping oscillation; 


as used herein, <£— tan 1 



phase angles between lift or pitching 
moment and displacement associ- 
ated with jth mode of oscillation, 
deg; for example, 


(pL. Qj = tan 1 


(CL.Q f 


angular frequency, radians/scc 
angular frequency in jlh mode, ra- 
dians/sec 


imaginary and real components 
mode of vibration under considera- 
tion 

chordwise and spanwisc pressure 
modes, respectively, in aerody- 
namic quantities (see oq. (13)) 
indicates parabolic camber deforma- 
tion 

indicates pitching 
indicates flapping 


THE INTEGRAL EQUATION RELATING LIFT AND 
DOWNWASH 

The linearized boundary-value problem of de- 
termining the forces on a wing with either oscil- 
latory or steady down wash conditions can be 
readily reduced to a problem of solving an integral 
equation that relates downwash and lift distribu- 
tion. The purpose of this section is to introduce 
this equation and to cast it in a desirable working 
form. 

PLAN FORM AND COORDINATE SYSTEM 

Since the integral equation is derived in various 
publications (e.g., in refs. 20 and 21), the equation 
will not he rederived here but will be simply stated 
in terms of dimensionless coordinates that arc 
particularly convenient for use in applications 
pertaining to flutter. In keeping with linear 
theory, the wing is considered as a plane impene- 
trable surface S that lies nearly in the xy-planc as 
indicated in sketch 1. In accordance with the 
common practice in analytical flutter studies of 
defining displacements and forces as positive 
downward, the ~-axis is taken to he positive 
downward. Furthermore, the x,y,z coordinate 
system and the surface S are assumed to move 
in the negative ^-direction at a uniform velocity 



V. It is remarked that the direction of the 
s-axis differs from that employed in oilier treat- 
ments (e.g., ref. 21) of the integral equation and 
its kernel function. This difference will be seen 
to lead to changes in sign in the kernel and in the 
relation between downwash and displacement 
when this procedure is compared with the treat- 
ments in references 14, 22, and 23. 


THE INTEGRAL EQUATION 

With the considerations mentioned in the pre- 
ceding section, a convenient form of the integral 
equation (obtained, e.g., from eq. (1) of ref. 21) 
may be formally written in terms of the unknown 
pressure distribution Ap(£,7j,t) as 

w(x,y,t) hi C l 
V 4xpV*J-i 

P' ^(Z,y,t)K[M,k,x 0 ,s(y--r))]d£dri (1) 

where w(x f y,f) is the vertical velocity or down- 
wash at any point (x,y) on the wing, p is the fluid 
density, and K[M f k } XQ } ft(y—ij)] is the kernel 
function of the integral equation. The quantity 
Kjiirp V (having the dimension of feet per second 
per pound) is the mathematical expression for the 
contribution to the downwash induced at any 
point (/, ?/) due to a unit force acting at any other 
point ({, if). In reference 21, K has been reduced 
to forms suitable for numerical evaluation. In 
equation (1), k is the reduced frequency b Q o>jV] x 0 
denotes the quantity x— £ in which x and £ are 
dimensionless chordwise variables referred to the 
root semichord b Q ; V and ?? are dimensionless span- 
wise variables referred to the semispan l y a con- 
venient reference length in dealing with the span- 
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wise integration (note that this differs from the 
usage of ref. 21 in which y and rj arc referred to 
b 0 ); s is the ratio of semispan to root semichord 
l/b Q ; and £ le (v) and £, c (??) are the chordwise co- 
ordinates of the leading and trailing edges, re- 
spectively. 

A solution to the integral equation for a given 
plan form, a known mode of oscillation, and a given 
set of stream conditions requires a determination 
of the pressure distribution which will satisfy the 
edge conditions appropriate to the plan form and 
flow regime under consideration and which, when 
multiplied by the kernel function and integrated 
over the plan form, will yield the down wash distri- 
bution corresponding to the mode of oscillation. 
No means of solving the integral equation in an 
exact analytic sense has yet been found and it has 
therefore been necessary to resort to approximate 
or numerical procedures. 

As stated in the introduction, a reliable concept 
on which to base an approximate solution to the 
integral equation involves expressing the unknown 
lift distribution as a sum of appropriately chosen 
modes of lift functions, each mode being weighted 
by a constant coefficient which must be determined. 
Employing this concept then allows the known 
down wash distribution to be expressed as a sum 
of definite integrals with the unknown coefficients 
appearing as factors of these integrals. Once the 
definite integrals arc evaluated, the unknown co- 
efficients can be readily determined by simple 
collocation. 

In order to cast equation (1) in a useful working 
form it is convenient to consider first some desired 
forms of the various ingredients of the equation, 
namely, the down wash distribution, the lift dis- 
tribution, and the kernel function. 

FORM OF DOWNWASH DISTRIBUTION 

In order to obtain a desired form of the down- 
wash distribution for use in equation (1), it is 
assumed that the system under consideration is 
undergoing a displacement TI(x,y,t) which may 
be represented by a superposition of either natural 
or assumed modes of vibration so that 

HU',y,t)=hi(r,y) qi(t)+h 2 (r,y)q 2 (t)+. . . 

+h J (r,y)q J (t)+. . . (2) 

where, for sinusoidal oscillations, 

(3) 


specifies the magnitude of the displacement in 
the jth mode, co is the frequency of oscillation, and 
hj(x,y) gives the shape of the mode. 

The dowmvash w(x,y,t) associated with the dis- 
placement II(x 7 y,t) is given by 

(4) 

or, with the use of equation (2), 

+h»(r,y) ^4- • •] (5) 

In application, each mode shape appearing in 
equation (5) will be handled separately. For con- 
venience the downwash associated with the jth 
mode may be defined as 

( 6 , 

FORM OF LIFT DISTRIBUTION 

The choice of the lift functions selected to ex- 
press the unknown pressure distribution should, 
of course, be made with discretion because the 
more compatible these functions arc with the 
actual loading the fewer w ill be required. As a 
guide in selecting modes of lift, the character of 
the lift distribution, at least in the neighborhood 
of the wing edges, can he surmised from a knowl- 
edge of the few exact solutions to lifting-surface 
problems. For subsonic flow, the pressure or lift 
distribution on a wing should tend to zero along 
trailing edges and edges parallel to the stream 

direction as lim ->[e } where € is the distance to the 
«->o 

wing edge. At the leading edge, the pressure or 

lift should behave as lim ^p. A summation of 
«->0 \ 6 

modes that satisfies these conditions and is otherwise 
perfectly general is employed, for example, in refer- 
ences 10, 13, and 14 and is used in the present 
development. 

Before introducing this summation of modes, 
it is convenient to introduce an angular chordwise 
variable 6 , shown in sketch 2. 

The functions £«(»?), and £ t e(v) in 

sketch 2 represent equations of the leading edge, 
midchord line, and trailing edge, respectively. It 


qj(t)=qje iut 
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Sketch 2. 


is noted that they may apply to plan forms with 
curved edges as well as to those with straight 
edges. The chordwise coordinate £ may be ex- 
pressed in terms of b/b 0 , and 0 through the 

relation 

£=£ m — -r-cos 0 (O^ 0 ^tt) (7) 

Note that and b/b d can be readily expressed in 
terms of %u(y) and as 


^rn 


£rg~l~ %le 
2 


b tu-Zic 

b 0 2 


( 8 ) 


The assumed pressure distribution of references 
10, 13, and 14 employs the angular chordwise 
variable 9 and is of the form 


&p(Z,y,t)=&p(O,v,t)=4*pV 2 frL(0,n)%jp- (9) 

VO U 0 

where 

T I— Q 

£(0,77) =-^Vl — y 2 \ (^00 ~h V a oi ~F V 2 ^02 T" • • d Cot g 

+ (a l0 +mn+v 2 ai 2 + . . .) sin 0+ . . . 

“b^n ( ff 7io+^ni - h , ? 2a »24" • • .)sin n$-\- . . .J 

( 10 ) 

in which the coefficients a nm are unknown weight- 


4 

ing factors to be determined. (The factor °f 

sin n6 , not emplojed in the cited references, is 
used here for convenience.) 

Equation (10) can he written more concisely 
by grouping the chordwise and span wise pressure 
modes to obtain 


in which 




h(B) = cot ^ (i?=0) 


and in which 


(ID 


(12) 


/!„(*?) = Vl — r) 2 ( a nO-\~V a ni J r « • • • • ♦) 

(13) 


Note that in dealing with symmetric loading, 
only even powers of 77 are retained in equation 
(13); in dealing with antisymmetric loading, only 
odd powers of t\ are retained. For either sym- 
metric or antisymmetric loading, as will be ob- 
served in the collocation procedure to be used, it 
is necessary to consider control points on only one 
wing panel. 

It is remarked further that equations (7), ( 8 ), 
and ( 12 ) may be used to obtain the following 
corresponding algebraic forms of the chordwise 
pressure modes: 


h(e)=m=j^l (£) 2 -(u -$) 2 

W)= 


1 ( 14 ) 


I) 2 


and so forth. For application to the special case 
of surfaces with straight leading and trailing edges, 
the following definitions of £* c and £, e iti terms of 
the span-to-chord ratio s, the taper ratio X, and the 
sweep angle of the quarter-chord line A may be 
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useful: 


— l + h/|(* tau A+-~) 

f «. = 1 + \v\ f an A — 3 ^ 2 ~ X -j 

Equation (0) and (lie associated definitions 
constitute a useful form of the assumed pressure 
distribution. In order to proceed toward a 
working form of equation (1), if remains to con- 
sider the form of the kernel function. 

FORM OF KERNEL FUNCTION 

A satisfactory working form of the kernel 
function is obtained from the expression for this 


function that is given by equation (D8) of refer- 
ence 2 1 . For later considerations it is convenient 
to have the kernel function expressed in a sep- 
arated form as 


/ffM,/-, Jo, «(?/—!?)] 


*«■**(?/— n>* 


05) 


where K denotes a modified form of the kernel 
function obtained by extracting from K t lie factor 

ft 2 s< 2 (y_ 7? ) 2 ' ma .V he noted, this factor can 

give rise to a second-order singularity and, hence, 
necessitate use of the concept of the “finite part 
of infinite integrals” in integrations involving the 
kernel function. For the oscillatory case the 


modified kernel function K is defin 


ed as 


A'f A f,k,.r t „x(i/-r))]- er 


ik*\y— 1)[ -) jj] K i (/ .s'| ■?/— nD + 'I ^I-V — v I [ h (£*!?/— i?l) — T‘\ (Hv— »)[)] + 


To- 


j r 5 * Y i h- v i 

-iks\y—ri] 

J 0 




Vi +e 


e < k *W-v\r ( l T 


(16) 


where Ki and /, arc modified Bessel functions and 
/-i is a modified Struve function. For the steady 
case (£=0) the modified kernel function K is 
defined as 

K[M, 0, r 0 , s(y jj)]=1 07) 

Examination of K shows that it contains no 
infinite singularities but does possess one finite 
singularity when x Q =y— r? = 0. This singularity 
leads to no special difficulty, however, since it can 
he handled^by making use of the following limiting 
forms of K: 


Win K[M ) k i T ii ,8(y—Ti)] = 2e- 1 * r * 

y — r )->0 

lim K[Af } k y ro,s(y—r))]=0 


Co>0) "I 


Co<0) 


y (is) 


Equation (1.5) and the associated definitions of 
equations (16), (17), and (18) constitute the 
desired working form of the kernel function for 


use in the integral equation. (Note that in the 
definition of the kernel the sign differs from that 
given in ref. 21. This is a consequence of the 
convention used herein that forces are positive 
downward.) 

Before returning to the integral equation it is 
desirable to discuss means of evaluating some of 
the functions which appear in equation (1C). 
First, the integral will be treated; then, the Bessel 
function K x {hs\y~ q|) and the combined function 
I\{ks\y— v\)~ L\ (frs|.y — i?l) will be discussed. 

The integral in K,— It lias not been possible to 
evaluate the integral which appears in equation 
(10) in an exact analytical sense. Although this 
integral is of a simple type for numerical integra- 
tion, its evaluation can become time consuming 
or burdensome, particularly when the upper limit 
becomes large. In the present procedure an 
approximate, but accurate, alternative approach 
lias been adopted. 

By examination of equation (10) it may he seen 
that if x o y>Afs\y—7j\, the upper limit of the 
integral is positive, but if x 0 <iMs\y— tj|, the upper 


530273 GO 
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limit is negative. For present purposes it is 
desirable to express the integral in a form sueli 
that the integration extends only over positive 
ranges. The integral is thus expressed as 

P " (,r„> M»\y-v\) (10) 

Jo VH t 2 


The combined (pianlity Ii—Li can be replaced 
by a definite integral, a convenient form of which 
is obtained from a general expression for I, — L, 
given in reference 24 (p. 425). For the special 
case v=i this definite integral is 

h (H?/— 'll) — (k*\y— ’ll) 


and 



e~ m ^~ nU dr 

Vl+r 2 


(s 0 <M«\>J-v\) (20) 


where 


2/,-.d?/— 

X 



7T 



e -ks ]l/-,!cos « ( J a 

!\ — T 2 e~ ks \ v ~* T dT 

( 22 ) 


px\y- - v\ 

and 

, i \f\*Q 2 '+P 2 * 7 (y--y ) 2 
d 

Willi these forms, the approach is to approximate 
a part of the integrand by the following expres- 
sion, which applies for all positive values of t\ 


Either of the integrals in this equation can he 
readily evaluated by numerical integration meth- 
ods to any desired degree of accuracy or, alter- 
natively, an approximation to the integrand can 
he used which is integrable in closed form. The 
latter procedure has been employed in the present 
program because it is more economical of comput- 
ing time and is believed to be sufficiently accurate. 
It has led to the following expression: 

Ii{Hy-y\)~-L { {ks\y-y\) 


T 


Vi+r 2 


1 - o . 1 o 1 e - 1 ° 329r - 0 . S99<r 1 • 4067r 

— 0.09480933e _2 ' 9 ° 7 sin ttt 


( 21 ) 


The expressions which result when equation (21) 
is substituted into equations (19) and (20) arc 
readily integrable in closed form. The results 
have been incorporated into the computing pro- 
gram of the present procedure and need not be 
written here. 

The approximation given by equation (21) lias 

the same value as .. T - at the two limits r— 0 
VH-t 2 

ami r= 00 . The maximum deviation of the approx- 
imation from the value of , - „ is about 0.24 

V'1+t 

percent in the vicinity of r=1.5. 

The Bessel and Struve functions. - For values 
of k$\y — j? 1>0 the modified Bessel function 
K x (ks\ y i?|) is evaluated by use of series expressions 
contained in a computing routine provided through 
an organization called SHARE made up of users 
of the IBM 704 and 709 electronic data process- 
ing machines. When ks\y — i?| = 0 the product 
£% - n \KMy~ y\) has a limiting value of unity. 


2ks\ JtL 

x 


ul /_ 

\1.3 


1 .0085 ksly-rj] 


.3410+ 1 .0050 Fs*(y-» j) 2 
+[^—0.8675 ks]y-r,\ 


0.4648-1 0.9159 ks]y- y\' 




3410+ k 2 s*{y— 17 ) 2 


e~ kg]y 


(23) 


The error of this approximation reaches a maxi- 
mum of about 0.4 percent in the vicinity of 
kf<\y-~ rj | — 4 and is less for both smaller and larger 
values of the argument. The overall effect of 
this error in the whole problem is minor. _ 

An idea of the general character of the kernel K 
can be had by examining the curves shown in 
figure 1 . In this figure the real and imaginary 
parts of K are plotted against kx Q for several con- 
stant values of k«\y—r) j. The particular values 
shown are for a Mach number of 0 . 8 . 1 be signifi- 

cant features of these curves are the finite discon- 
tinuity at k*\y r}\ — kxQ~0 and the rapid variations 
in K near A\r 0 = 0 especially for small values of 
£% 77 ]. As will be seen, this behavior of the 

kernel necessitates special handling of the chord - 
wise integrations over a range of values of ks\y- y\ 
near y rj—O. 
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(a) Real part. 

Figure 1 . — Variation of modified kernel function i?|) with kx Q 

for various values of 1 ?|. M =0.8. 


WORKING FORM OF THE INTEGRAL EQUATION 

The desired working form of the integral equa- 
tion (eq. (1)) is readily obtained by employing the 
various ingredients in the forms arrived at in the 
preceding sections. By making use of the down- 
wash as expressed by equation (5), the lift dis- 
tribution as expressed by equation (9), and the 
kernel function as expressed by equation (15), the 
integral equation may be written for thcjth mode 
of oscillation as 

sin 6 dO (24) 

where 

K 9 (M i k J B,ri)=K^M i k, [x—tm+y cos ??)J 

As was remarked earlier, the presence of the 


second-order singularity necessitates the 

use of the concept of the finite part of infinite inte- 
grals indicated in equation (24) by the notation 

The concept, as here employed, is discussed, for 
example, in references 11 and 25. Thus, the value 

of the improper integral f for a<?/<6 is 

Jn ( y v ) 

obtained from the following limit: 


7 r F(v)dv 

.(y-O* ™Uo iy-vY 

+£.£pSH fW ] < 2 ® 


Equation (24) together with equation (25) con- 
stitutes the desired working form of the integral 
equation. The subscripts j on the modal function 
h and on the lift function L in equation (24) indi- 
cate that a solution for the pressure distribution 
must be obtained for eacli individual mode of 
oscillation. In the discussions to follow, this will 
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(b) Imaginary pari. 
Figure 1. Concluded. 


be understood to be the case and the subscripts 
will be dropped. 


SOLUTION OF THE INTEGRAL EQUATION 

When the expression for L(6,r}) given by equa- 
tion (10) or (1 1) is substituted into equation (24), 
the downwash is expressed as the following sum of 
definite integrals: 


w 


V n m J -I t, 


/l — n 2 <h 


j; 


(y-v) 2 

l n (d)Ke (MJkfitf) sin Odd (26) 


The next problem, then, is one of devising a prac- 
tical scheme for accurately evaluating the inte- 
grals. 

The scheme of integration which has been de- 
vised is simple in principle but because of the 
many operations and parameters involved, it rep- 
resents such a task of calculation that it is feasi- 
ble only with the use of high-speed computing 
equipment. It is well suited for such equipment, 


however, and once it is programed the approximate 
lift distribution or its various integrated properties 
can be generated for a series of given mode shapes 
in a very few minutes. 

In discussing the procedure if is convenient to 
separate the chordwise and spamvise integral 
operations and to write equation (2G) as 

w(r,y)_ C l F(yj)ibe) (?7) 

V J-i (y—y ) 2 

where 

F(v) s W T^/nW (28) 

n 77i 

The function /«(i?) in equation (28) is the chord- 
wise integral and is defined by 

f n (v) - In ( d)Ke(M,k , e, v ) Sin e <le (29) 

EVALUATION OF THE CHORD WISE INTEGRAL 

The integrals defined by equation (29) involve 
no infinite singularities and can therefore be readily 
evaluated by numerical means to any desired de- 
gree of accuracy. In selecting specific integration 
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procedures to be employed, the clmracteristics of 
the integrand should, of course, he kept in mind 
and provision should he made for handling certain 
special situations. Specifically, proper account 
must be taken of the finite discontinuity in the 
kernel function at ./' 0 =2/ — ^7=0, indicated by equa- 
tion (18), and of the rapid variation in the kernel 
near x 0 = 0 for small values of 1 ?|, indicated 

in figure 1. 

In the present procedure these characteristics 
of the kernel are taken into account by performing 
the chord wise integration in two steps for certain 
ranges of values of k«\y~r)\. For this purpose 
equation (29) is written as 

fn(v)= j /„(0)i^(A/,/i*,0,?7) sin 6 dO 

+ | l n {6)Ke{M i kfi ) ri) sin OdO (-10) 

where 0 X denotes the value of 0 which corresponds 
to x 0 — 0. The value of 0 Z to he used in equation 
(30) can he obtained from equation (7) as 

^os-R (*"-')] (-'5 1 ) 

By closely spacing the integration stations in 
each integral of equation (30) in the vicinity of 
the end point 0 X (as is done, e.g., in a Gaussian 
procedure), the rapid variation in the kernel 
can be more accurately handled. 

An appropriate range of values of jfc*|y— rj\ over 
which to employ the two-part chord wise integral 
(eq. (30)) can he arrived at by examining plots 
like those of figure 1. In the present procedure 
equation (30) is being handled by two 10-point 
Gaussian integrations in the region bs\y— 0 . 3 . 
For larger values of ks | y — rj\ } the chord wise 
integrals are being evaluated by use of equation 
(29), with a single 10-point Gaussian integration. 
The computing program of the present procedure 
is extremely flexible with regard to the number 
of integration stations employed and can he easily 
modified. 

One remark concerning the application of 
equation (30) to swept configurations is perhaps 
pertinent. When a control point is near a leading 
or trailing edge, it is possible for the line %=x 
or Jq—O to intersect the edge within the range of 
7) values in which the two-part chord wise integra- 
tion is used. In this event the argument of the 


inverse cosine in equation (31) will exceed ±1.0. 
In the computing program of the present procedure 
provision is made for testing the argument of the 
inverse cosine and for performing the chord wise 
integration in one part, by use of equation (29), 

wherever ~ \^ m — ^[>1.0. 

Some sample values of chord wise integrals 
obtained by the present procedure are shown in 
figure 2. The figure shows the products of the 
factor \ \—rf and the chord wise integrals / {J (??), 



FinuRio 2. — Variation of ehordwise integrals f n (y) (eq. 29) 
will) y for a control point at y~ 0.5 on a swept tapered 
wing. ,17 = 0.81; k - 0.3; A = 3.3. 
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f\(v) 7 and /ofa), defined by equations (29) and 
(30). Tliese products are shown as functions of y 
for a tapered wing with straight edges. Note 
that for this ease there is a discontinuity in slope 
of the curves at y~ 0 which is due to the factor 
bjb in the lift function L(0,y). Tliis factor is the 
ratio of root chord to local chord and, for the 
type of plan form treated, is discontinuous in 
slope at the wing root . 

It is remarked that in earlier unpublished 
versions of the present procedure the chord wise 
integrals have been expressed in terms of the 
algebraic form of the chordwise pressure modes 
7 n (£), defined by equation (14). The first mode 
/^(£) contains a singularity at the leading edge 
which can be handled by subtraction and addition 
in the integrand of the corresponding chordwise 
integral. This is a straightforward step but it 
leaves the integrand with infinite slope at the 
point where the singularity existed. As is well 
known, for integration of such shapes, closely 
spaced integration stations in the vicinity of the 
infinite slope arc needed to approach exact results. 
Use of the angular variable 0 , together witli 
Gaussian integration procedures, should yield 
more accurate results for a given number of 
integration stations. 

THE SPAN WISE INTEGRATION 

In order to accomplish the spanwisc integration 
it is convenient to divide the wing into several 
regions as indicated in sketch 3, 



Region I extends from 17 = — 1 to 77=0 and is 
employed to take into account a possible discon- 
tinuity in slope at 77=0 which may arise from the 
factor b 0 /b as indicated in the preceding section. 
Region III extends a short distance f on each side 
of y, that is, from y = y— f to 17 — y+f- The 
integrand of equation (27) for this region contains 


the second-order singularity ^ __^ 2 ; hence, this 

region is intended to facilitate the evaluation of 
the finite part of the improper integral. Region II 
simply fills the gap between regions T and III and 
is absorbed into region III in cases for which the 
control point is in the neighborhood of 77 = 0 . 
Region IV fills the gap between region III and 
77 = 1,0 and is absorbed into region ITT in cases for 
which the control point is near a wing tip. Note 
that the function F(y) in regions I and IV is 
characterized by an infinite slope at the wing tips 
which arises from the factor Vl —y 2 . 

In conformity with the four regions just dis- 
cussed, the spanwisc integral can be written as 


where 


f" f(n)>lv 

' ' J , (v-vY 

_ p t F( qV/q 

( y-vY 

r v+t F(ri)di) 
-r ( y-vY 

j f 1 F(y)dy 
1Y ~)y + X (1/ VY 


(3311) 


(:3:5b) 


(33c) 


(33d) 


First, the evaluation of the integrals I i, In, and 
/ IV will he discussed. Then, the integral 7m, 
which contains the singularity, will be treated. 

The integrals /,, 7 n , and 7 IV are not singular 
and can be readily evaluated by numerical means. 
In the present, procedure the function F(rf) at 
each integration station within a particular region 
is evaluated from equation (28), divided by the 
appropriate value of (y—yY, and weighted by the 
integrating factor for the station. Thus, the 
integral I 1} for example, is expressed as 


.=i 


F(v>) 

iy—y*Y 


(34) 


where the constants A, are the integrating factors 
and indicates the coordinates of the integrating 
stations appropriate to the Appoint integration 
rule being used. 

With Tegard to appropriate values of N, if is 
remarked that in the present procedure the span- 


AERODYNAMIC FORCES ON OSCILLATING OR STEADY FINITE WINGS AT SITBSONIC SPEEDS 13 


wise integrals are also handled by Gaussian inte- 
gration rules. 4 lie number of integration stations 
employed varies from region to region and depends 
on the spamvise position of the control points and 
on aspect ratio. For aspect ratios of the order of 4 
or less, region I is being handled by a 10-point rule 
regardless of control -point location; region IT is 
being handled by from a 3- to a 10-point rule, the 
number increasing as the control point is moved 
outboard; in region IV, even though the control 
point maybe well outboard, at least a 5-point rule is 
used to handle the aforementioned infinite slope at 
the tip and the number of points is increased to 10 
as the control point moves inboard. For wings of 
much higher aspect ratio the number of points used 
in all three of the regions should be increased. 
As in the case of the chord wise integrals, the span- 
wise-integration scheme of the present, procedure 
is vcrj r flexible with regard to the number of 
integration steps employed and the number can 
be easily changed. 

The integral 7 IU , which involves the singularity 
at 7) = y } is evaluated by employing a polynomial 
to approximate F(rj) over the region of integration. 
This procedure permits the singular parts of 7 m to 
be completely isolated and easily handled. In the 
present scheme F(tj) is represented by a sixth- 
degree polynomial that is readily determined by 
use of Lagrange’s interpolation formula and is 


S P 0? 2/+35)(?y— ?/+2<$)(q— y-\ 6) 

W -p)l V~y+(3-p)i 5 

X (n — ay) (v—y—&)(y—y—26) (v-y— 36) F p + X 

(35) 

where 5 is the station interval (<5=f/3) and F p +i 
denotes the value of the chord wise integral at sta- 
tion ^+1. These quantities are indicated in 
sketch 4. 



integral /, u can be written as 





F(v)</ r} 

(y-vY 


1 rv+t 

+6%(,V-V)+& 4 <U} </v+\ Ur, ' r — 

<5 J v-iy-v 

f»+t (In 

(30) 

where each term g is a certain linear combination of 
the ordinates of F(rf) and can be found by use of 
equation (35). For the present discussion only 
the final combination of the ordinates of F(rj) is of 
importance and the expressions for g need not be 
written. 


The notation (j^ in equation (30) is meant to 

imply that a Cauchy principal value of an inte- 
gral is to be taken. For the integral involved the 
Cauchy principal value is 

*L + r™ *L\ 

J v-ty—v €^o\Jv-3sy— n Jv+t y—yJ 

=lim[Iog 35— log «— log (—35) 

e ->0 

+log(— c)]=0 (37) 

TJie last integral in equation (3G) involves the 
singularity and is to be evaluated by use of the 
concept of the finite part of infinite integrals. 
Tie nee, 

f +f jL =lim r c— <h | r +n <h n 

J v-t(y vY c ->0 L^/i/-35(2/ — ^) 2 Jv + t (y— 1 ?) Z ej 

_ii nl (_!_|- + _!_r +, ‘_?V = 2 

U- 35 y—y \v+t e/ 3<$ 

(38) 

Tlie nonsingular integral in equation (30) yields 

i r v+ * 

? J f f \{h^y~Tty j r6(ji(y~riY 


+8 2 g 2 (y-vY+^(y~v)+s i g 4 ](h 


U =y~- (729<7 o + 1 3% 2 + 45;/.) (39) 

Tn order to accomplish the integration nse is 

made of the identity y=y—(y—y) so that t lie Thus, combining equations (37), (38), and (39) 
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gives the integral Tm fis 

I m «^^(729</(H- 13503+450*— 5f/ a ) (40) 

This result may be expressed in terms of the 
ordina tes P\ t . . - Pi of 7* ( 17 ) HR 

-^ s (. 3 F 1+ 72F, + 4 9.^- 1 ,300 F< 1 49 r>F 3 

+727^ 6 +l 3/^) (411 

The complete spamvise integral is given by the 
sum of the results indicated by equation (34) for 
T u T u , and 7 1V and given by equation (41) for I lu . 
This sum expresses the dowmvash at tlie control 
point under consideration as a linear algebraic 
equation with the weighting factors a„ m as un- 
knowns. 

in order to determine the values of the unknown 
weighting factors a nm , the simplest procedure is to 
select just as many control points as there are un- 
knowns and to form and solve the set of simul- 
taneous equations related to these points. 1 his is 
an area where further work could he done in deter- 
mining whether there is an optimum number of 
control points for a given problem or whether 
something like a least-squares approach might he 
useful. These questions are not settled here; some 
further remarks on the need for numerical research 
arc given in the next section. 

REMARKS ON THE METHOD 

The purpose of this section is to indicate briefly 
some areas in winch the need for further exploia- 
tory calculations and developments exists and to 
provide some remarks on the use of the present 
program. 

areas for further numerical research 
and development 

"Phe development of the method presented 
herein has left untreated some questions and 
problems which may arise in applications and 
considerations of the procedure. Among these 
are the questions of optimum numbers or locations 
of control points or of the number of pressure 
modes necessary to assure satisfactory results for 
a given case. 

Unfortunately, it is not possible, at least at 
this time, to give firm answers of general appli- 
cability to such questions. It appears unlikely 
that a general optimization of either control-point 
location or number of pressure modes is feasible 


hut that each plan form and dowmvash distribution 
must bo considered as a separate problem. It is 
hoped that for general applications optimization 
will not be necessary hut that satisfactory accura- 
cies can he obtained by employing a relatively 
few pressure modes along with a well-scattered 
pattern of control points. 

A great deal seems to hinge on the accuracy with 
which the chosen lift functions represent the actual 
lift distribution for a particular problem. As the 
functions chosen to represent the lift become more 
and more realistic, the locations of control points 
should become less significant. The lift functions 
employed in the present procedure arc chosen with 
a hope that the forms employed are at least near 
enough to reality to allow control points to he 
scattered in a pattern which is adequate for 
representing the dowmvash distributions en- 
countered. Applications of the procedure to the 
circular (or elliptical) plan form, for which some 
results of analytic methods exist, might serve as 
the basis for studies of such questions. 

While t lie computing program as it now exists 
is directly applicable to controls such as all- 
movable tails, there is a need for modifications 
and extensions that will make it applicable to 
controls in general. This would involve the 
determination of appropriate lift functions for 
handling various edge conditions which may he 
involved in control problems. 

More experience is needed in applying the 
procedure to problems involving higher modes of 
oscillation and the associated higher frequencies. 
Presumably, for such problems, greater numbers of 
control points and lift functions would ho needed. 
This, of course, implies an increase in the order of 
(lie matrices to he handled and may lead to for- 
mation of less well conditioned or even poorly 
conditioned matrices. 

Some study is also needed to determine the 
benefits, if any, of (lie application of methods of 
weighting, such as least squares, that might 
permit evaluation of the unknown coefficients 
from a large number of control points rather than 
from the minimum number required. 

COMPUTING PROCEDURE 

As mentioned previously the procedure described 
heroin lias been programed for the IBM 704 
electronic data processing machine. As presently 
programed, provision has been made for the 
inclusion of either three or four ehordwise lift 
functions, together with either three or four 
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spanwisc lift functions, in tlic series given by cqua- 
tion (10). This programing permits direct use of 
cither 9, 12, or 16 control points. 

The program proceeds in two steps: In one step 
the surface integrals associated with each control 
point and, thus, the coefficient matrix of the 
unknowns in the problem are generated; in the 
second step the coefficient matrix is inverted and 
combined with a number of sets of specified 
dowmvash conditions to obtain the weighting 
fac tors a nm associated with each downwash con- 
dition. The two steps are separate so that should 
it be desired to apply, for example*, a least-squares 
approach, as many equations as desired can be 
generated and an appropriate means of solving 
them provided. 

In most of the calculations performed thus far, 
the usual practice has been to include three 
cliordwise terms and three spanwisc terms of 
equation (10) and to use these with just 9 control 
points. Usually the control points have been 
placed at or near the >(-, ){-, and %-chorcl stations 
and at approximately the 0.2-, 0.5-, and 0.8- 
semispan stations. In (he calculations described 
herein for the circular plan form, a very limited 
number of control-point variations have been 
attempted. Calculations have been made with 9 
control points located as just described and with 
9 control points arrayed as recommended in 
refc'i ence 1 1 . Calculations have also been made 
with 12 control points by adding 1 control point 
at each of the three span stations used in the 9- 
point calculations and by distributing the points 
cliordwise in the manner recommended in refer- 
ence 11. Tn this limited study, results obtained 
f°i J dl coefficients varied by less than 2 percent. 

The calculations consume about 25 seconds of 
machine time per control point when the numeri- 
cal integrations are carried out with about the 
number of integration steps indicated in the text. 
Thus, about 4 minutes are required for the calcu- 
lations with 9 control points. This time includes 
both program steps described in this section and 
involves the application of about five sets of 
downwash conditions in the second step. 

The basic information provided by the comput- 
ing program is a set of values of the pressure 
weighting factors a nm for each prescribed down- 
wash distribution. The program can be readily 
extended by the introduction of appropriate inte- 
grating matrices, and with very little increase in 
machine tune, to yield any desired integrated 
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property of the pressure distribution such as sec- 
tion or total lift and moment coefficients or the 
generalized forces required for flutter analysis. 

NUMERICAL APPLICATIONS 
In order to illustrate applications of the method, 
results of several calculations arc presented. 
(Some of the results are taken from ref. 2G; other 
results are as yet unpublished.) In the illustra- 
tions total forces and moments are compared (1) 
with results of analytic procedures for a circular 
plan form with steady downwash conditions, (2) 
with results of other theories and with experi- 
ments for a rectangular plan form of aspect ratio 1 
at a uniform angle of attack, and (3) with some 
experimental results for a rectangular plan form 
of aspect ratio 2 undergoing pitching and flapping 
oscillations. Also included arc results of flutter 
calculations compared with experiments for an all- 
movable control surface of aspect ratio 3.50 and 
for a cantilevered rectangular plan form of aspect 
ratio 5.04. 

THE CIRCULAR PLAN FORM 

Of the lew plan forms for which analytical 
methods of treatment have been developed, one is 
the circular plan form. A compilation of results 
for the total lift and pitching moment on a circular 
plan form in steady incompressible flow" is given 
in reference 9. The present procedure has been 
applied to this configuration with tw'o distribu- 
tions of downwash. These are for (1) an unde- 
formed surface at constant angle of attack, de- 
fined by h(x,y) =x, and (2) a surface with parabolic 

camber deformation, defined by A(.r, y)=~ 

Total lift and pitching-moment coefficients and 
the location of the center of pressure have been 
obtained by the method presented herein and com- 
pared with results listed in reference 9 from the 
various analytic procedures. For the undeformed 
surface at constant angle of attack the results are 
as follow's: 


Source 

C t .a 

Cm, a 

Center of 
pressure, 
percent 
chord 
from lead- 
ing edge 

Khmer (ref. 27) _ 
Schade and Kricnes 

L 8174 

0. 9358 

24. 75 

(ref. 7) 

1. 7984 

. 9318 

24. 00 

Van Spiegel (ref. 9)__ 

1. 7902 

. 9320 

23. 95 

Present method 

1. 7910 

. 9389 

23. 80 
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For the surface with parabolic camber deforma- 
tion the following results have been obtained : 


Source 

C L ,J 

Cm, 

Center of 
pressure, 
percent 
chord 
from lead- 
ing edge 

Tvinner (ref. 27) — 
Sehade and 

TCrienes (ref. 7)__ 
Van Spiegel (ref. 9). 
Present method 

0. 9350 

. 9436 
. 9326 
. 9443 

-0. 4376 

4382 

4388 

4463 

20. 00 

20. 80 
20. 45 
20. 40 


The results of the present method differ by less 
than 2 percent from the results of the analytic 
procedures. All the results of the analytic pro- 
cedures have been obtained by summing the first 
few terms of infinite series and, hence, are not 
exact. Were more terms of the series taken, or 
more significant figures carried in the calculations, 
the last figure shown (and possibly the last two 
figures) of the results of the analytic procedures 
might change. 


AERONAUTICS AND SPACE ADMINISTRATION 

THE RECTANGULAR PLAN FORM OF ASPECT RATIO 1 

In reference 28 results from various theories 
have been collected and compared with an extrap- 
olated experimental result for the total lift and 
center of pressure for a rectangular wing of aspect 
ratio 1.0 at a constant angle of attack. One of 
the tabulated results, obtained by use of the initial 
computing program of the present procedure, 
agreed poorly with the experimental result. In 
the initial program chordwise integrations were 
performed by use of Simpson’s rule with rather 
widely spaced intervals. In addition, the chord- 
wise lift distribution was expressed in algebraic 
form and the leading-edge singularity was ex- 
tracted analytically. Such usage generally results 
in a small loss of accuracy arising from the infinite 
slope of the integrand at the leading and trailing 
edges. In this case there seems to have been a 
sensitivity to the loss of accuracy. The current 
computing program employs improved integration 
techniques and yields a considerably better result, 
as shown by the following comparison with some 


of the results of the various procedures listed in reference 28 


Source 

! 

| 

Center of 
pressure, 
percent 
chord 

Experiment (extrapolated) (ref. ^ 

Jones (zero-aspect ratio) (ref. ~9) -- 

1. 400 
1. 570 
1. 400 
1. 497 
1. 310 
1. 455 

16. 00 
0 

17. 00 
14. 37 
8. 79 
17. 21 

Lawrence (ref . 30) _ 

ii.su trei. - - — — "* 

Present method (initial program) 

Present method (current program) 


THE RECTANGULAR PLAN FORM OF ASPECT 
RATIO 2— COMPARISON OF THEORETICAL 
AND EXPERIMENTAL AIR FORCES 

Tlie results of a number of experiments on 
simulated two-dimensional wing models (wings 
with end plates) arc in good agreement with 
theoretical two-dimensional oscillating air forces, 
hut relatively few experimental programs have 
been carried out to measure the air forces on oscil- 
lating finite wings. Among the few which do exist 
is a systematic series of tests on rectangular semi- 
span wing models of aspect ratio 2. Two modes 
of oscillation have been studied; one, a pitching 
oscillation about an axis along the midchord, and 


the other, a flapping oscillation about an axis 
inboard of the wing root and parallel to the stream 
direction. The tests were conducted by two 
coordinated techniques with semirigid spring- 
mounted models. One technique was to force the 
model to oscillate; then, the resonant frequency 
was determined, the power input to the oscillator 
was obtained, and the lift forces were measured 
with strain gages. These values were compared 
with tare values determined from a resonant oscil- 
lation under vacuum conditions. The other 
technique was to oscillate the model at its resonant 
frequency, cut off the applied force, and record 
the decay of the motion. By testing over a range 
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of airspeed, results were obtained over a range 
of reduced frequency. 

The pitching wing. — The experimental results 
for a midchord location of the pitch axis, obtained 
in reference 31, arc presented in figure 3 and are 
compared with theoretical results for the kernel- 
function procedure. Figure 3(a) gives the ampli- 
tude and phase angle of the total lift coefficient 
C Lf a plotted against the reciprocal of k, the re- 
duced frequency, and includes a sketch showing 
the semispan model, the pitch axis, and the tunnel 
wall. Figure 3(b) gives the amplitude and phase 
angle of the total pitching-moment coefficient 
Cm, a- Also given along the abscissas of figures 
3(a) and (b) are the test Mach numbers corre- 
sponding to various values of 1 JJc. The theoretical 
curves were obtained by using the appropriate 
Mach number. The agreement of theory and ex- 
periment is generally good but the results in figure 
3(b) show appreciable scatter in the data. 

The flapping wing. — In figure 4 the results from 
reference 22 for a flapping mode of oscillation are 
shown. Amplitudes and phase angles for the 
total lift coefficient C LA and the total pitching- 
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(a) Total lift coefficient. 


Figure 3* Comparison of theoretically and experi- 
mentally determined lift and pitching-moment coeffi- 
cients (amplitudes and phase angles) plotted against 
1/fc for a rectangular wing of aspect ratio 2 oscillating 
in pitch about midcliord. 


STEADY FINITE WINGS AT SUBSONIC SPEEDS 17 



\/k 

(b) Total pitching-moment coefficient.. 
Figure 3. — Concluded. 


moment coefficient Cm# (midchord axis) appear 
in figures 4(a) and] 4(b), respectively, plotted 
against the reduced frequency k. Figure 4(a) also 
includes a sketch of the semispan model, tunnel 
wall, and the flapping axis behind the tunnel wall. 
The test Macli number ranged only up to 0.41 and, 
therefore, the theoretical curves were determined 
on the basis of A/== 0, The agreement of the- 
oretical and experimental coefficient amplitudes 
is good and the phase-angle agreement ranges from 
good to fair. A considerable decrease in scatter 
of the experimental data in comparison with the 
data of figure 3 is evident and is believed to be 
the result of improved experimental technique 
and instrumentation achieved over a period of 
time. 

COMPARISON OF THEORETICAL AND EXPERIMENTAL 
FLUTTER RESULTS 

A number of applications of the kernel function 
procedure have been made in calculating flutter 
boundaries, but relatively few published compari- 
sons with experiment have been made available. 
One such comparison appears in reference 23 for a 
45° delta semispan model mounted as a cantilever 
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(a) Total lift coefficient. 

Figure: 4. — Comparison of theoretically and oxperi- 
men tally determined lift and pitching-moment coeffi- 
cients (amplitudes and phase angles) plotted against k 
for a rectangular plan form of aspect ratio 2 undergoing 
a flapping oscillation. 


O Experiment, reference 22 
Theory, M - 0 



at its root and tested at M— 0.85. In this study 
the theoretical boundary was found to be slightly 
conservative compared with the experimental 
boundary. Another comparison of theory and 
experiment, indicated in reference 26, is based on 
a scries of flutter experiments carried out in the 
Langley 8-foot transonic pressure tunnel. The 
model, shown schematically" at the top of figure 5, 
was a swept, tapered, all-movable control surface 
with an aspect ratio of 3.50. In figure 5 the 

flutter-speed ratio t — ^—j= and frequency ratio w/co 2 

are plotted against M. The model mass ratio m 
equals the mass of the model divided by the mass 
of the air in a truncated circular cone with the 
local diameter equal to the local chord of the model. 
The mass ratio n is included in the flutter-speed 
ratio because wind-tunnel density had to be grad- 
ually changed in order to obtain flutter at different 
Mach numbers on the same model. The second 




k 


(b) Total pitching-moment coefficient. 
Figure 4. — Concluded. 



igure 5— Theoretically determined flutter-speed ratio 
^ - and flutter-frequency ratio «/co 2 plotted against 

0|)0>2y M 

M and compared with two flutter experiments on an all- 
movable control-surface model. (See ref. 26.) A = 20 , 
AO: \ = 0.15. 



AERODYNAMIC FORCES ON OSCILLATING OR 

natural frequency u 2 is used because the second 
natural mode most nearly resembled a pure pitch- 
ing mode. The theoretical results shown were 
obtained by using three nat ural modes. The speed 
ratio is unconservative by a few percent in com- 
parison with the two experimental points shown; 
the flutter frequencies agree quite well. 

Another comparison between theoretically and 
experimentally determined flutter characteristics 
is shown in figure G for a cantilever-mounted rec- 
tangular wing of aspect ratio 5.04. The experi- 
mental data in this figure represent a portion of 
(he results of a series of tests presented in reference 
32. In figure G the stiffness-altitude parameter 
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— and frequency ratio afa are plot led against 

M. For the configuration of this example the 
second natural mode of frequency is essentially 
a pure torsion mode. The theoretical results 
shown were obtained by using the first three un- 
coupled modes of a uniform cantilever beam. For 
the range of Mach number over which comparisons 
are made, the theoretical values of the stiffness- 
alt. tude parameter agree very well with the ex- 
perimental values; the theoretical values of the 
flutter-frequency rat io are within 10 percent of the 
experimental values, 

concluding remarks 



Fif iURE 0. -Theoretically determined id iff n ess-altitude 
parameter --VJI and flutter-frequency ratio 

plotted aganist M and compared with flutter experi- 
ments on a cantilever-mounted rectangular plan form of 
aspect ratio A of 5.0 1. 


A detailed description has been given of the 
numerical procedure employed in a solution of 
the mteg m l equation which relates oscillatory or 
steady lift and downwash distributions on finite 
wings m subsonic flow. The procedure, which 
has been programed for the IBM 704 electronic 
data processing machine, is applicable to general 
plan forms with either curved or straight leadin'- 
and trading edges. While, as developed, the pm! 
ccdure is readily applicable to control surfaces 
such as all-movable (ails, modifications are needed 
to apply it to controls in general. 

The procedure has been applied to a number of 
examples, some of which deal with the determina- 
tion of integrated aerodynamic forces and moments 
and some of which deal with applications to flutter 
problems. Integrated forces and moments ob- 
tained for a circular plan form are compared with 
results of analytical procedures for some steady 
downwash conditions; the results differ by less 
than 2 percent. Integrated forces and moments 
lor rectangular plan forms are compared with ex- 
periment or with other theories for steady flow 
and for pitching and flapping oscillations. In each 
instance the results are in generally good am-ec- 
mail. ' * 

The flutter examples include comparisons of 
calculations with experiments for an all-movable 
control surface and for a rectangular wing of aspect 
ratio 5.04, both at high subsonic speeds. For both 
examples the theoretical results agree well with 
tile experimental results. 


Langley Research Center, 

National Aeronautics and Space 
Langley Field, Va., June 2, 1959 . 
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